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Abstract 

We study the fractional populations in chromosome inherited diseases. The 
governing equations for the fractional populations are found and solved in the 
presence of mutation and selection. The physical fixed points obtained are 
used to discuss the cases of color blindness and hemophilia. 

87.10.+e, 42.62.Be 



1 



Typeset using REVTgX 



I. INTRODUCTION 



Many physical ideas are finding applications in complex biological systems these days 
Recently, we presented a theoretical scheme in which one can investigate the ratios 
between the fractional population of blood groups ||. This method has some analogy with 
the physical concept of renormalization and fixed points. In this paper, we extend the theory 
to investigate the problem of sex-linked inheritance. 

In sex-linked inheritance, there are five population groups, XX, X'X, X'X', XY, and 
X'Y, where X, X' and Y represent the normal female, the defective female, and the male 
chromosome, respectively. The group of so called carrier female is characterized by X'X. 
The genetic rule is that sons receive Y and daughters do X or X' from their fathers. Similarly, 
their mothers deliver X or X' to sons and daughters. It is also known that there exist 
mutations between X and X' . Another important factor in the present problem is that the 
defective groups, X'X' and X'Y, have disadvantages in surviving and inheriting unlike in 
the case of the blood groups [[J . This selection process should be taken into account for any 
reasonable discussions. Therefore, we consider the inheritance of sex-linked disease in the 
presence of mutation and selection and obtain the governing equations, which determine the 
next fractional populations from the previous ones. The governing equations will be used 
to investigate the problems of genetic propagations of chromosome-linked diseases such as 
color blindness and hemophilia. 



II. FRACTIONAL POPULATION EQUATIONS 

We consider the five fractional populations with the following constraints: 

XX(n) + X'X(n) + X'X'(n) = 1 (1) 
XY(n)+X'Y(n) = 1, 

for the n-th generation. The ratios of gene frequencies without mutation can be determined 

as 

Xf(n) = XX(n) + -X'X(n), 



X' f (n) = X'X' {n) + -X'X(n), (2) 



1 
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X^(n)=XY(n), 
XL(n) = X'Y(n). 



Although the mutation between the normal chromosome X and the defective one X' 
would be rare, still it plays an important role in the following discussion. In order to 
consider the mutation, we introduce two probability factors, a and j3 for the following 
mutation processes; 

X X' and XJ-X'. (3) 
Through the above mutation processes, the gene frequencies are modified as 
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X f (n) = (1 - a)X f (n) + /3X' f {n), 

X' f (n) = (l-(3)X' f (n) + aX f (n), (4) 

X m (n) = (1 - a)X m (n) + /^(n), 
I» = (l-/3)X;( n ) + «I m (n). 

The fractional population equations, which govern the populations of the next generation, 
are now written as 

XX{n + l)=X f {n) xX m {n), 
^(n + 1) = X' f (n) x X m (n) + X f (n) x X^(n), 
A^'(n + l)=Xj(n)x^(n), (5) 
XY(n + l)=X f (n), 
X J Y{n + 1) = Xj(n). 

Since the defective chromosome causes a disease, the populations of X'X' and X'Y will 
have less chances of surviving and inheriting their genes. In order to reflect this disad- 
vantages, we introduce disadvantage factors Sf for the female and 5 m for the male groups, 
respectively. Then, the populations of X'X' and X'Y will be modified as 

PoplX^X'in + 1)] — ► (1 - 6 f )Pop[X 7 X'(n + 1)], (6) 
Poplin + 1)] — ► (1 - 6 m )Pop[JCY(n + 1)]. 



With normalization, the fractional population equations are given by 

XX(n + l) = 



X f (n) ■ X m (n) 



X'Xin + 1} 



l-5 f X' f {n)-X' m {ny 
X' f (n) ■ X m (n) + X f (n) ■ Xjjn) 
l-5 f X' f {n)-X' m {n) 



, . (1 - SAX'Jn) ■ X'(n) 



XY(n + l) 
X'Y(n + 1) 



Xf(n) 



l-6 m X' f ( n y 
(l-5 m )X' f (n) 
l-5 m X' f (n) 



The above governing equations yield the following constraint relations for any gen- 

eration n, 

X f (n)+X' f (n) = l and X m (n) + X' m (n) = 1. (8) 

In order to understand the change of populations along generations, it is convenient to 
consider the automata equations for X'j- and X' m only, which are given by 
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mi + ] ~ 1 - 5 m X' f (n) ' (9) 



2a + (1 - a - /?) (x}(n) + X^(n)) - 2fr(l - /?)X}(n) • X^(n) 
/( " +) " l- W (n).X^(n) 

The above coupled recursion relations can now be used to study the fixed points of X' m 
and XL which correspond to the equilibrium values where X' m (n + 1) = X' m {n) = X' m * 
and X'j(n + 1) = X'j(n) = X'^* . The fixed points of X'f* are given by the solutions of the 
algebraic equation, 

a + ai X' f * + a 2 {X' f y + a 3 (X' f *) 3 = 0, (10) 
where the coefficients are given by 
a = a (3 — a — ft), 

a x = (a + ft){a + ft - 3) + 5 m (2p - a - (3 2 - a(3 - 1) + 25 f a{(3 - 1), (11) 
a 2 = Ml + " + P) + 25/(2a + 2(3 - (3 2 - a(3 - 1) + 25 m 5 f ((3 - l) 2 , 
a 3 = 25/(1 - a -p) + 26 m 5 f (p - 1). 

Solving this equation for stable fixed points, we can readily determine the equilibrium pop- 
ulation ratios. 

We study the fixed points in several cases. First of all, in the case of no mutation; 
a = and ft — 0, the meaningful fixed point, X'j*, is given by 0. It correctly predicts 
that without mutations, the defective genes will disappear eventually. Secondly, when only 
mutations are considered in the theory assuming 5 m = and 5f = 0, the fixed point is given 
by X'j* = a I [a + ft). This result will be used in the discussion for color blindness. In other 
general cases, the exact solution cannot be expressed in a closed form. However, since the 
mutation rates are known to be very small(~ 10~ 5 ~ 10~ 7 ), the fixed point X'^ can be 
expressed in terms of a and ft in an approximate fashion and will be discussed in the next 
section. 



III. COLOR BLINDNESS AND HEMOPHILIA 

The disadvantage, that a color blindness man or woman has, is not severe enough to 
reduce the chance of survival significantly. Hence, we let the disadvantage factors be simply 
zero in this case. Then, we easily notice from Eq. (|9p that the fixed point is given by 

This result is identical to that obtained in the conventional genetics Q . 

Using the fractional population equations of Eq. (0) and the above fixed point, we obtain 
the population ratios as 

XY : X'Y = X) : X' f * = 1 - X'/ : X' f * = ft : a, (13) 
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and furthermore 



XX : X'X : X'X' = (3 2 : 2/fa : a 2 . (14) 

The above result is the well known the Hardy- Weinberg law 0. 

The demographic data for color blindness in England show that XY : X'Y = 12 : 1 ||. 
Hence we conclude that /3 : a — 12 : 1. We notice that X' chromosome is more unstable 
than X chromosome since (3 is much larger than a. The abundance of carrier female is easily 
noticed by the fact that XX : X'X : X'X' = 144 : 24 : 1. 

Hemophilia is a dreadful disease which affects the chance of survival and mating signifi- 
cantly. All of the female group X'X' perish completely upon birth. Hence, the disadvantage 
factor 5f is equal to 1. Then the fixed point X'j* can be expressed up to the second order of 
a and (3 as follows using Eq. fllOD and ([TT|), 



8 m \ 8 m S m 8^ J \8i 

3 )a+ ( 2 10 + 30 15 

8 m \ 8m 8™ 8™ I \ 8m $L 



^* = (-2 + ? )«+f-2-^ + ^-^V 2 +f-2 + ^-ll^ + 0(a 3 ; ...). (16) 



"m m/ \ m m/ 



A numerical calculation of the stable fixed point, X'** using Eq. ([H]) is shown in Fig. 1. We 
find that the dominant contributions come from a and 8 m as Eq. (^) indicates. We also find 
that the values of the fixed point Xj* are independent of initial values. The approximate 
expressions, Eqs. ([15]) and (0) are found in good agreement with the exact expression Eq. 
(prO|) except when 8 m is near zero. 

It is useful to consider the male hemophilia population before selection in order to relate 

the above formulation with the statistical data. Here, the relevant demographic data is the 

— * - — - * 

ratio, r, between the mutation cases and the all hemophilia cases; X'Y mut = r ■ X'Y . Using 
Eqs. (@)~(0), the population of male infants with hemophilia before selection can be written 

as 

XV* = \(1 + « - P)X'X* + aXX* = JCY* nh + XV* mut . (17) 

Here, the second term represents the male population having hemophilia caused by the 
spontaneous gene mutation. Actually, the first term also contains the gene mutation contri- 
bution of \{a — (3)X'X* . However, when the demographic data r are collected, there is no 
way to distinguish the inheritance from the mutation in this case, because data collectors 
simply check whether there were hemophiliac occurrences in the family line or not. It is 
straightforward to show that r is related to the disadvantage factor 8 m and the mutation 
rates a and j3. From Eqs. @~(0), we find 

5 m = 3r + (10 - lOr - -)a + (-3 + 5r)/3 + 0(a 2 , a/3). (18) 

The early statistical data for the infant male population having hemophilia show that 
XY : X'Y ~ 10 4 : 1. Assuming that the current fractional population distribution has 
reached a fixed point, we find X'j* ~ 10 -4 from Eq. (||). Furthermore, a recent statistics 
shows that the rate r is given by r ~ | ||. This data and Eq. yield 8 m = 1 — 0(a, (3). 
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This result is in a reasonable agreement with the fact that various therapies treating male 
hemophilia have been invented only recently and, thus, the probability of successful marriage 
and reproduction was almost zero for male patient in the past. With 6 m ~ 1 and X'j* ~ 10~ 4 , 
we find from Eq. fllED that the value of a is about 3.3 x 10~ 5 , which is in a reasonable range 
as the probability of mutation. Also, we can obtain the population ratio of the carrier female 
group using Eq. (|;XI: X'X : X'X' ~ 10 4 : 1.33 : 0. 

Since the level of therapies treating male hemophilia has now reached the stage that most 
of male patient can marry and reproduce, it is interesting to study the case 5 m = and make 
predictions how the mutation rate and the fractional population will change accordingly. The 
fixed point of Eq. flTSP and ([TBI) will be modified as 
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Using the above results, we find the modified rate 



X? = Jl a $ + 0(a,-.-) t (19) 



X' m * = a + (1 - a + (3)X' f \ (20) 



r= P a h+0(a,---). (21) 

Assuming the mutation rate a does not change and remain a ~ 3.3 x 10~ 5 as we find in the 
above, we can determine X'^* = 0.00707. Also, we readily obtain the fractional populations: 
XX : X'X : X'X' = 98.6 : 1.41 : 0, and XY : X'Y = 993 : 7. The result clearly 
shows that the majority of hemophilia would results from inheritance and that the fractional 
population of carrier female increases drastically, when male hemophiliac patients survive 
and mate without any disadvantages. Therefore, nongenetic treatment of hemophiliac male 
may cause increase of hemophiliac population and infant deaths of female patients, unless 
some concurrent measures are taken. However, we note that the present calculation can not 
produce the dynamic properties of the transition period between 5 m ~ 1 and S m ~ 0, since 
5m is assumed static in the calculation. 



IV. CONCLUSION 

We have considered the population ratios of the genetic groups related with chromosome 
inherited diseases. The governing equations, which determine the ratios, are found in the 
presence of mutation and selection. The selection is taken into account in the formulation 
by using the disadvantage factors. It is found that there exist physical fixed points in the 
automata equations, which correspond to equilibrium population rates. These fractional 
population equations are used to discuss the cases of color blindness and hemophilia. 

In the case of color blindness, there is no significant disadvantage in selection so that the 
disadvantage factors can be assumed zero. From the governing equations, we readily obtain 
the Hardy- Weinberg relation XX : X'X : X'X' = (3 2 : 2a(3 : a 2 . Using the statistical data 
XY : X'Y = 12 : 1, we find that the ratio of mutation rates a : /3 = 1 : 12. 

Hemophilia seriously hampers chances of survival, mating and reproduction. Especially 
for female patients, chance of survival is almost zero, thus, making 5f = 1. Using the 
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demographic data that one out of ten thousand male infants have hemophilia, and that one 
third of all hemophiliac cases are thought to be caused by gene mutation, we obtain the 
following results; i) the disadvantage factor 5 m for male is almost 1, ii) the mutation rate 
a ~ 3.3 x KT 5 , and hi) XX : X'X ~ 10 4 : 1.33. 

We have also studied the case when the hemophiliac male suffers no disadvantage in the 
selection process; 5 m = 0. It is found that the population of hemophiliac females and males 
would increase drastically and inheritary hemophilia would be dominant over gene mutated 
cases. 
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FIGURES 



Figure 1 




FIG. 1. For a set of values of 5f = 0.5 and (3 = 10~ 5 , we plot the fixed point X'j* in the 
three dimensional format, where the x and y-axis correspond the mutation rate a and the male 
disadvantage factor 5 m . It is found that the overall features of the shape and size do not depend 
on 8f and (5 sensitively. 
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